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Abstract. The light curve of the quasar OJ 287 extends from 1891 up today without ma- 
jor gaps. This is partly due to extensive studies of historical plate archives by Rene Hudec 
and associates, partly due to several observing campaigns in recent times. Here we sum- 
marize the results of the 2005 - 2010 observing campaign in which several hundred scien- 
tists and amateur astronomers took part. The main results are the following: (1) The 2005 
October optical outburst came at the expected time, thus confirming the general relativis- 
tic precession in the binary black hole system. This result disproved the model of a single 
black hole system with accretion disk oscillations, as well as several toy models of binaries 
without relativistic precession. In the latter models the main outburst would have been a 
year later. (2) The nature of the radiation of the 2005 October outburst was expected to be 
bremsstrahlung from hot gas at the temperature of 3 x lO'' °K. This was confirmed by com- 
bined ground based and ultraviolet observations using the XMM-Newton X-ray telescope. 
(3) A secondary outburst of the same nature was expected at 2007 September 13. Within the 
accuracy of observations (about 6 hours), it started at the correct time. Thus the prediction 
was accurate at the same level as the prediction of the return of Halley's comet in 1986. (4) 
Further synchrotron outbursts were expected following the two bremsstrahlung outbursts. 
They came as scheduled between 2007 October and 2009 December. (5) Due to the effect 
of the secondary on the overall direction of the jet, the parsec scale jet was expected to 
rotate in the sky by a large angle around 2009. This rotation may have been seen at high 
frequency radio observations. OJ 287 binary black hole system is currently our best labora- 
tory for testing theories of gravitation. Using OJ 287, the correctness of General Relativity 
has now been demonstrated up to the third Post-Newtonian order, at higher order than has 
been possible using the binary pulsars. 
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1. Introduction 
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In 1982 Aimo Sillanpaa put together the his- 
torical light curve of OJ 287 based on pub- 
lished measurements. There appeared to be 
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Fig. 1. The optical light curve of OJ 287 from 
1891 to 2010. It includes unpublished data 
from R.Hudec and M.Basta 



a 12 year outburst cycle (see Figure 1), and 
moreover, it was obvious that the next cyclic 
outburst was due very shortly. Indeed, OJ 287 
produced the expected event in the following 
January (Haarala et al. 1983; SiUanpaa et al. 
1985). Observations showed a sharp decline 
in the percentage polarization during the out- 
burst maximum, indicating that the outburst 
was produced essentially by unpolarized light 
(Smith et al. 1987). This is different from or- 
dinary outbursts in OJ 287 which are charac- 
terized by an increase of the percentage po- 
larization at the maximum light. Radio out- 
bursts were found to follow the optical out- 
bursts with a time delay of between 2 months 
and a year, depending on the observing fre- 
quency (Valtaoja et al. 1987). 

Sillanpaa et al. (1988) suggested that 
OJ 287 is a binary black hole system where 
a smaller companion perturbs periodically the 
accretion disk of a massive primary black hole. 
The next expected outburst was in 1994; it 
came as scheduled (?Sillanpaa et al. 1996a). 
In the binary model there should be two disk 
crossings per 12 yr orbital period. Thus the 
1994 outburst should have an equal pair whose 
timing was calculated to be at the beginning 
of 1995 October (Valtonen 1996; Lehto & 
Valtonen 1996; SundeUus et al. 1996). This 
prediction was also verified (Sillanpaa et al. 
1999b). 

Alternative explanations were also put for- 
ward. Quasiperiodic oscillations in an accre- 



tion disk were suggested (Igumenshchev & 
Abramowicz 1999), and several binary toy 
models without relativistic precession were 
also proposed (Katz 1997; Villata et al. 1998; 
Valtaoja et al. 2000). The latter models all pre- 
dicted the next main outburst of OJ 287 in the 
autumn of 2006 while the precessing binary 
model gave a prediction one year earlier, at 
the beginning of 2005 October (Sundelius et 
al. 1997; Kidger 2000). The second major out- 
burst was expected in late 2007 in all binary 
models while in the oscillating accretion disk 
model there was no reason to expect a second 
major outburst. In the precessing binary model 
the date was given with high accuracy, with the 
last prediction prior to the actual event being 
2007 September 13 (Valtonen 2007, 2008a). 
In the oscillating accretion disk model and in 
non-precessing binary models the nature of 
the radiation at these outbursts should have 
been polarized synchrotron radiation while the 
precessing binary model predicted unpolarized 
bremsstrahlung radiation (Lehto & Valtonen 
1996). In addition, the precessing binary model 
predicted a series of further outbursts for the 
interval 2007 - 2010 (Sundelius et al. 1997). 
In this model the companion black hole should 
also affect the disk of the primary in a pre- 
dictable way, leading to the wobble of the jet 
(Valtonen et al. 2006). 

With these predictions in mind, a multi- 
wavelength campaign of observing OJ 287 dur- 
ing 2005 - 2010 was set up. 

2. Five 'smoking gun' results 

In the following we describe five 'smoking 
gun' observations which produced expected re- 
sults from the point of view of the precessing 
binary black hole model but which were sur- 
prising and unexpected in other theories. 

2. 1. Timing tlie 2005 outburst 

The 2005 outburst was well covered by alto- 
gether 2329 observations in V and R-bands. 
The points in Figure 2 are daily averages, 92 
in all. 

According to (Sundelius et al. 1996), the 
impact causing the 2005 outburst was expected 
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Fig. 2. The optical light curve of OJ 287 during 
the 2005 outburst. The dashed line is the theo- 
retical fit based on (Lehto & Valtonen 1996). 



22.3 years after the impact of the 1983 out- 
burst. In addition, (Lehto & Valtonen 1996) es- 
timated that the 2005 outburst should be de- 
layed after the impact. Also the 1983 outburst 
is delayed but not as much, the difference be- 
ing 0.44 yr. The rapid flux rise started in the 
latter outburst at 1983.00; thus the correspond- 
ing rapid flux rise of the 2005 outburst was ex- 
pected at 1983.00 + 22.30 -I- 0.44 = 2005.74. 
Actually, the outburst was one week late and 
did not begin until 2005.76 (Valtonen et al. 
2008c), but anyway the timing was well within 
the error limits. Only a few polarization mea- 
surements were carried out at that time, and 
unfortunately, even those happened during sec- 
ondary flares. Thus the polarization state of the 
primary outburst is not known. 

The timing confirmed the precession of 
39.1° per period. It is so much higher that e.g. 
in binary pulsars (by a factor of 10"*) that we 
immediately realise the importance of OJ 287 
in testing General Relativity. 

2.2. Nature of radiation at ttie 2005 
outburst 

The impact outbursts are expected to consist 
of bremsstrahlung radiation and thus the op- 
tical polarization of OJ 287 should go down 
during them. As mentioned above, the polar- 
ization information for the basic 2005 outburst 
is not available. However, bremsstrahlung may 
be recognized also by it spectrum, and this is 



Fig. 3. The optical - U V spectrum of OJ 287 
during the 2005 outburst. The soUd Une is the 
bremsstrahlung fit as predicted in (Lehto & 
Valtonen 1996). The observational points are 
corrected for the internal extinction in OJ 287 
as well as for the extinction in our Galaxy. 



w polarization flux in OJ287 




Fig. 4. The optical light curve of OJ 287 during 
2006-2008. Only low polarization (less than 
10%) data points are shown. 



the part of the campaign that was successfully 
carried out. 

We had XMM-Newton observations both 
before the 2005 outburst (2005 April), and 
during the outburst (2005 November 3-4). 
Fortunately the November observation hap- 
pened at the time when the source was at its 
basic outburst level, in between two secondary 
bursts. Thus we would expect to see an addi- 
tional pure bremsstrahlung spectrum above the 
underlying synchrotron power-law. A prelimi- 
nary report of these observations has appeared 
in (Ciprini et al. 2007), and a more detailed re- 
port is under preparation. 
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Fig. 5. The optical light curve of OJ 287 during 
the 2007 September outburst. Only low polar- 
ization (less than 10%) data points are shown. 
Data points are based on Ref. 22. The dashed 
line is the theoretical fit based on Ref 9. The 
arrow points to September 13.0, the predicted 
time of origin of the rapid flux rise. 



below 10% (Valtonen et al. 2008b). Thus it is 
not difficult to decide which one was the ex- 
pected bremsstrahlung event. Later in the year 
there were more highly polarized outbursts, 
but if we look at the Ught curve composed 
of low polarization states only (Figure 4), the 
September 13 outburst clearly stands out as a 
sharp spike. The light curve has been plotted 
with a better resolution in Figure 5. A theo- 
retical light curve is also drawn, and an arrow 
points to the expected moment of the beginning 
of the sharp flux rise. We see that the observed 
flux rise coincides within 6 hours with the ex- 
pected time. The accuracy is about the same 
with which we were able to predict the return 
of Halley's comet in 1986! 



2.4. 2007 -2010 outbursts 



In Figure 3 we show the difference between 
the 2005 November flux and the 2005 April 
flux. The values have been corrected for the 
Galactic extinction and for the extinction in the 
OJ 287 host galaxy. For the latter, we use the 
measurements in ? and the standard Galactic 
extinction curve (Draine 2003). The soUd line 
shows the bremsstrahlung spectrum at the ex- 
pected temperature of 3 x 10^ °K. Note that 
a raised synchrotron spectrum, as one might 
have expected in some other theories, would 
have a downward slope toward higher frequen- 
cies, and it is entirely inconsistent with obser- 
vations. 

2.3. Timing and nature of tiie 2007 
September 13 outburst 

The 2007 September 13 outburst was an ob- 
servational challenge, as the source was visible 
only for a short period of time in the morning 
sky just before the sunrise. Therefore a coor- 
dinated effort was made starting with observa- 
tions in Japan, then moving to China, and fi- 
nally to central and western Europe. 

Four outbursts of various sizes were de- 
tected in 2007 September. In three of them the 
degree of polarization was above 15% while 
the fourth and the biggest one had polarization 



Sundelius et al. (1997) gave a detailed predic- 
tion of the whole light curve of OJ 287 dur- 
ing the campaign period. In addition to the 
two impact outbursts, it was expected that the 
tidal forcing mechanism of (Sillanpaa et al. 
1988) would create a long hump in the fight 
curve, starting in the middle 2007 and contin- 
uing until the middle of 2009. This hump is 
clearly seen in Figure 5, together with the 2007 
September 13 spike. Seta et al. (2009) carried 
out multifrequency observations of OJ 287 at 
the hump (2007 November 7-9) and compared 
the flux values with earlier observations before 
the hump (2007 April). They find that the ex- 
cess radiation at the hump is neither due to in- 
creased magnetic field strength nor increased 
Lorentz factor, but likely results from an in- 
crease in electron energy density. This is what 
is expected in the tidal model. 

A prominent outburst was also expected at 
the end of 2009 and it was observed (Valtonen 
etal. 2011a). 



2.5. Jet reorientation 

The accretion disk as a whole is also affected 
by the companion in its 12 year orbit. If the 
jet and the disk are connected the jet direction 
should be influenced by the companion. 



286 



Valtonen: OJ 287 



Model 3 
+3yr 

dots: observations 


iv-'' . !• j 






\ 




I r\! 


1900 1950 2000 2050 



Fig. 6. The observed position angle of the radio 
jet of OJ 287 (points) compared with the binary 
model, with a 3 year jet response time to the 
reorientation of the disk. 



There are three periodicities that could be 
expected to show up: the 12 yr orbital cycle, 
the 120 yr precession cycle (or half of it due 
to symmetry) and the Kozai cycle (Innanen et 
al. 1997) which also happens to be ~ 120 yr. 
The 12 yr orbital cycle produces the tidal en- 
hancements in accretion flow, and the enhance- 
ment can be stronger or weaker depending on 
where we are in the precession cycle. These 
two tidal eff'ects explain the overall appearance 
of the light curve (Sundelius et al. 1997). In 
addition, there is a modulation in the long term 
base emission level due to the jet wobble which 
follows the Kozai cycle. 

The jet wobble shows up in observations in 
several ways. First, the mean angle of optical 
polarization varies (SiUanpaS et al. 1991). The 
binary model predicts, among other things, a 
quick change in the optical polarization angle 
by nearly 90° around 1995 which has been ob- 
served (Villforth, C. et al. 2010). In radio, we 
should see a similar rapid change in the posi- 
tion angle of the parsec scale jet. Depending 
on the delay in radio jet reorientation, the 
change could be already under way (Figure 6), 
or it may be delayed by another 12 yr cycle 
(Valtonen et al. 2011b). There are recent ob- 
servations which suggest the first alternative 
(Agudo et al. 2011), but the interpretation of 
these observations is not yet clear-cut. 



3. Testing General Relativity 

Using the OJ 287 binary, we may test the idea 
that the central body is actually a black hole. 
One of the most important characteristics of a 
black hole is that it must satisfy the so called 
no-hair theorem (Misner et al. 1973). A prac- 
tical test was suggested by Thorne (1980). In 
this test the quadrupole moment Q of the spin- 
ning body is measured. If the spin of the body 
is S and its mass is M, we determine the value 
of q in 



For black holes q = 1, for neutron stars and 
other possible bosonic structures q > 2 (Wex 
& Kopeikin 1999; Will 2008) 

We have studied the distribution of q- 
values in orbits which give the correct timing 
of 9 outbursts within the range of measurement 
accuracy. We find that the distribution peaks at 
q = I, thus confirming the no-hair theorem. 
The parameter values q = or q = 2 can 
be rejected at the 3 standard deviation level at 
present. This is the first time that it has been 
possible to study general relativity at higher 
than the 1.5PN order. 

4. Conclusions 

Prior to the 2005 - 2010 multiwavelength cam- 
paign there were several ideas about the na- 
ture of OJ 287. Fortunately, these models made 
completely different predictions about the be- 
haviour of OJ 287. The result of a scientific 
enquiry is seldom as clear-cut as this: the out- 
bursts satisfied the expectations of the precess- 
ing binary model in every respect while the al- 
ternative ideas failed. Therefore we are confi- 
dent that OJ 287 can be used as a test labora- 
tory for theories of gravitation. Our first results 
strongly support general relativity as the cor- 
rect theory. 

Acknowledgements. We would like to thank all the 
participants of this campaign for an extraordinary 
amount of work dedicated to the solving of the rid- 
dle of OJ 287. 
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DISCUSSION 

FRANK RIEGER: (1) What is the gravita- 
tional lifetime of the system? I presume it is 
short. (2) You seem to need that the secondary 
crosses the primary disk twice per orbit. 1 re- 
member that (Ivanov et al. 1999) found that the 
co-aUgnment (orbital plane vs. disk plane) hap- 
pens relatively quickly. Did you check for that? 

MAURI VALTONEN: (1) 10000 yr (2) 
Observational evidence is that the disk and the 
orbital plane are at roughly 90 degree angle rel- 
ative to each other. This comes primarily from 
the so called "eclipses" which may arise when 
the secondary crosses the jet. Theoretically, the 
first order perturbation theory shows that the 
incUnation can remain at 90 degrees if it has 
obtained a value close to it at some earlier time. 
Our numerical simulations confirm this. 

MANEL PERUCHO: There is tendency to 
associate the radio jet helical structures to pre- 
cession caused by binary black hole systems. 
Could comment on this point? Is the periodic- 
ity you find in the jet of OJ 287 the expected 
one from your model? 

MAURI VALTONEN: Marscher & Jorstad 

(2011) have discovered that in the large scale 
X-ray and radio jet there are oscillations in the 
time scale of thousands of years. In our model 
the spin axis of the primary precesses with the 
period of about 1300 years (Valtonen et al. 
20010). These two facts may be cormected. 



